Supermassive Black Holes from Primordial Black Hole Seeds 
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The observational evidence for a population of quasars powered by supermassive black holes of 
mass > 10 9 Mq at redshifts z > 6 poses a great challenge for any model describing the formation 
of galaxies. Assuming uninterrupted accretion at the Eddington limit, seed black holes of at least 
lOOOM© are needed at z~ 15. Here I study whether seeds could be primordial black holes (PBHs) 
which have been produced in the very early universe by the collapse of primordial density fluctu- 
ations. In particular, I study the expected number densities of PBHs in the relevant mass range 
for several classes of spectra of primordial density fluctuations and confront the results with obser- 
vational data. While it seems to be possible to produce the required PBHs with spectra showing 
large enhancements of fluctuations on a certain scale, our hypothesis can be clearly disproved for a 
scale free spectrum of primordial fluctuations described by a power-law slope consistent with recent 
observations. 

PACS numbers: 98.80 Bp, 98.62 Lv, 04.70 Bw 



I. INTRODUCTION 

Supermassive black holes (SMBHs) are strongly be- 
lieved to dwell in the centers of most galaxies. They are 
also thought to be the central engines of active galactic 
nuclei (AGN) and quasars |l|. While the above state- 
ments seem to be common knowledge, by now, the pro- 
cess of formation and subsequent evolution of these ob- 
jects is one of the fundamental problems of contemporary 
astrophysics. 

The problem concerning the seed black holes, which 
eventually grow to SMBHs, has sharpened during the 
last few years by some both accurate and interesting ob- 
servations. First of all, in the local universe a strong 
correlation has been found between the mass of a SMBH 
and the velocity dispersion of the bulge of its host galaxy 
• This provided evidence for the proposal that SMBHs 
have formed before and evolved together with the bulges 
of their host galaxies [3| and allowed to deduce comoving 
number densities for SMBHs of a given mass in the local 
universe. 

On the other hand, a population of quasars with red- 
shifts of z > 5.7 (this value of z marks the onset of the 
so called i-dropout) has been found by the Sloan Digi- 
tal Sky Survey (SDSS) 0, U @. Up to now twelve of 
these objects have been discovered and the most distant 
of them has a redshift z = 6.4 and a central black hole 
of mass M rs 3 x 10 9 M Q Q. This observation poses a 
great challenge for any model which seeks to describe the 
evolution of SMBHs from seed black holes by a combina- 
tion of mergers of dark matter halos, which host SMBHs, 
and accretion of baryonic matter. The merging history 
of dark matter halos can be described in the hierarchical 
model of structure formation according to the so called 
extended Press-Schechter formalism 71. However, after 



'Electronic address: nd@thp.uni-koeln.de 



a coalescence the newly formed SMBH gets a kick with 
some velocity which usually exceeds the escape velocity 
in a shallow dark matter halo. Thus the SMBH is cut 
off from its baryonic fuel supply sitting at the center of 
mass of the dark matter halo, and further growth of the 
black hole by accretion is no longer possible 8] . 

However, all models of growing SMBHs have to start 
off from a very early {z > 15) population of seed black 
holes, each of them having a mass of at least 1000 M©. 
Astrophysical suggestions for the production of these 
black holes comprise the collapse of the entire baryonic 
cloud at the core of a dark matter halo, the remnants 
of an early generation of very heavy and metal poor 
stars (Population III), or the formation of an Interme- 
diate Mass Black Hole (IMBH) from the collapse of an 
early star cluster. All these mechanisms which may be 
able to produce the seeds of the quasar population at 
z « 2 have serious problems to account for the black 
holes needed to have quasars as early as z rs 6 (For a 
recent review of this complex, see, e.g., |jj). 

As long as the origin of the mentioned black hole seeds 
is so unclear, one should have an open mind to alter- 
native scenarios. Here I study whether these seeds can 
be primordial black holes (PBHs) produced during the 
very first of cosmic evolution from the gravitational col- 
lapse of a sufficiently large overdense region 10] . For 
a given spectrum of primordial fluctuations, which may 
arise, e.g. , from some inflationary scenario, the expected 
mass spectrum of PBHs can be computed along the lines 
of, e.g., ref. [TjJ. The fluctuation spectrum is known to 
some degree for scales relevant for CMB and LSS obser- 
vations [l3, uM > but has to be extrapolated to smaller 
scales in order to compute the abundances of PBHs. It 
is the goal of this work to study whether the seeds of 
the z f» 6 quasars can be PBHs which arose upon the 
collapse of suitably parametrized primordial fluctuations 
consistent with recent observations. 

The paper is organized as follows: in the next section I 
describe shortly the problem of growing a SMBH by ac- 



cretion from a seed black hole. Then I shall discuss the 
formation of PBHs through the collapse of an overdense 
region in the very early universe, together with some pos- 
sible parametrizations of the spectrum of primordial cur- 
vature perturbations. In the following section, I study 
under which circumstances the resulting PBHs could be 
the seeds for the early population of quasars mentioned 
above. In the last section, I will discuss the results and 
draw my conclusions. 



II. ACCRETIONAL GROWTH OF SMBHS 

In this section, I give a simple estimation of the time 
needed to build a SMBH from a black hole seed by accre- 
tion. I will not go into details of the accretion process, 
but one should have in mind that the scenario presented 
here is highly optimistic, and so the calculated evolution 
times should be seen as lower limits to the real ones. 

Let us assume that a seed black hole of initial mass Mi 
sits at the center of some dark matter halo and starts at 
time ti to accrete radiatively baryonic matter with a rate 
Wfuci at the Eddington limit with an efficiency e. The 
luminosity of the radiation is then given by|33j 



L = em iucl c 2 = 1.15 -M. 

0~T 



(1) 



where m p denotes the proton mass, <jt the Thomson 
cross section for an electron, M, the (evolving) mass of 
the black hole, and the factor 1.15 is the mean atomic 
weight per electron for the hydrogen and helium gas mix- 
ture. Under the assumption that all accreted mass en- 
ergy, which is not radiated away, feeds the black hole|3J|, 
the SMBH grows according to the equation 



M, = (1 — e)rhf, 



ucl 



(2) 



which, upon inserting eq. (I), can be easily integrated to 
yield the solution 



t-U 



M.(t) = Mjexp 

where the e-folding time r is given by 
t = 1.24 x 



(3) 



xl0 16 s«3.9x — xl0 7 yrs. (4) 
1 : 0.1 " V ; 



The e-folding time is, of course, an increasing function 
of e, which itself is an increasing function of the angular 
momentum of the black hole and can reach values up to 
0.4 for a maximal rotating hole. So the value e « 0.1 of 
this so called thin-disc accretion mode, which is adopted 
here and also in a number of other papers on this sub- 
ject, seems to be quite optimistic for our goal to obtain 
short growing times for SMBHs. The assumption that 
thin disc accretion plays a major role in SMBH growth 
obtained further confidence recently by evaluating a rela- 
tion which is known as Soltan's argument 14]: integrat- 
ing the (measured) quasar luminosity function over all 



wave lengths and redshifts under the assumption of thin- 
disc accretion yields an estimate of today's mass density 
in SMBHs. This density turned out to be consistent with 
the density of the local SMBH population [lj| . 

Other accretion modes like the low effectivity advec- 
tion dominated accretion or the low efficiency super Ed- 
dington mode seem to suffer from considerably high out- 
flows, possibly connected to jets, or instabilities which 
prevents them from being effective for a sufficiently long 
time. However, super Eddington accretion may be im- 
portant for some part of the early accretion history |16J] . 
But even if this should be the case, because of large out- 
flows the process of SMBH growth is not expected to 
proceed faster than in the thin-disc mode. 

Our goal is to build up a SMBH of mass M, = 
3 x 10 9 Af Q at z = 6.4 from a black hole seed of mass 
Mi. As initial time for the accretion process we take 
a redshift z — 15, which is motivated by the onset 
of reionization 17]. For a flat ACDM-cosmology with 
Ho = 70km/(s • Mpc) and CIm — 0.3, we obtain a grow- 
ing time of t(z = 6.4)- t(z = 15) = 0.84 Gyr- 0.27 Gyr = 
0.57 Gyr. Then from eq. (J3J), we get a minimal seed black 
hole mass M, « 1300M Q . Here we have assumed that the 
accretion process is effective without any interruptions. 
It is often assumed that accretion occurs during certain 
cycles triggered by a larger merger (see, e.g., |l8|). These 
periods are then characterized by some duty cycle (usu- 
ally 10 7 — 10 8 years) or the amount of baryonic matter 
available for SMBH growth. Again our model seems to 
be quite optimistic, but one should also have in mind 
that at very early times the supply of baryonic matter is 
larger than after the onset of significant star formation. 

As mentioned above it is very challenging to provide an 
astrophysical mechanism to produce black holes of this 
mass at such an early time in the cosmic evolution |g|. 
The collapse of the entire baryonic cloud of a dark mat- 
ter halo seems to be problematic due to fragmentation 
after H 2 -cooling. So the first collapsed baryonic objects 
seem to be a population of superheavy (up to 1000 M Q ) 
stars, which, depending on their mass, collapse further 
to a black hole or explode via a pair instability super- 
nova without leaving any remants |l9| . The latter out- 
come would pollute the environment with metals, which 
trigger the more effective atomic line cooling mechanism 
and so prevent the formation of further superheavy stars. 
Even if the circumstances would favor the formation of 
black holes it seems to be very problematic to have them 
with the desired mass at such an early time (270 million 
years) of the cosmic evolution. 

This problem gives us the motivation for investigating 
whether the seed black holes for the formation of SMBHs 
could be primordial. 



III. PRIMORDIAL BLACK HOLES 

In this section we provide the formulae to calculate 
the mass spectrum for primordial black holes (PBHs), 



which are formed by the collapse of primordial density 
fluctuations, and discuss some ways to parametrize the 
spectrum of these fluctuations. 



PBH formation 



be safely taken to be infinity or, more convenient for nu- 
merical calculations, some value where the integrand is 
sufficiently small. 

Now we deduce the relations between various quan- 
tities for PBHs of mass Mjj(tk)- The horizon mass is 
defined to be 



Primordial black holes can be formed by the collapse of 
primordial density fluctuations |ifj. Because the Jeans 
radius in the radiation dominated era is of the order of 
the horizon scale, the mass of such a PBH is about the 
horizon mass Mh at the time of formation. Now the 
collapse of an overdense region is only possible, if the 
mis of the primordial fluctuations there, averaged over a 
Hubble volume, is larger than a threshold <5 m m- 

The value of 5 m ; n is and has been a matter of some dis- 
cussion. For a long time it had been thought to be about 
I /3. At the end of the nineties numerical simulations 
suggested that <5 min sa 0.7 would be more appropriate 
[20| . Very recently, an analytic calculation which employs 
peaks theory instead of modified Press-Schechter theory 
has been proposed which arrives at the result that both 
analytic approaches are in good agreement if one takes 
the Press-Schechter value <5 m in ~ 0.3 — 0.5 |21|. In the 
present work we take <5 m i n = 0.6 while keeping in mind 
that our results and constraints could be tightened or 
loosened by a different choice. 

A fluctuation mode with some comoving wave number 
k enters the horizon at time tk defined by k = a(tk)H(tk)- 
The probability for the formation of a PBH of mass 
Mfi(tk) at time tk is then given by 



P{M H ) 
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where ajj(t k ) denotes the variance of primordial density 
fluctuations at time t k , avaraged over one Hubble volume 
at that time. It is given by 



o%(tk) 
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dxx 3 A 2 R (kx)T 2 {kx,t k )W% H (x) 



(6) 

Here, A|j(fc) is the spectrum of primordial curvature fluc- 
tuations, Wth{x) is the Fourier transform of the top 
head window function given by 



Wth (x) = — r (sin x — x cos x) 



(7) 



and T 2 (k,t) is the transfer function for the subhorizon 
evolution of the density fluctuations. In the above inte- 
gral we are interested in the function T 2 {kx,t k ), which 
turns out to be W^. H (c s x), where c s — l/\/3 denotes the 
speed of sound in the radiation dominated era. The scale 
fee is an ultraviolet cut off for small scales which is needed 
if one chooses a fluctuation spectrum for which the inte- 
gral I© does not converge. Here one usually takes the 
scale of reheating. However, in the calculations of the 
present work the integral always converges, and k e can 
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(8) 

We are interested in PBHs of about 1000M©. From eq. 
(jHJ we see that these are produced well before the era of 
nucleosynthesis, but when we are dealing with fluctua- 
tion spectra we must ensure that no PBHs much heavier 
than the ones mentioned are produced, because the fluc- 
tuations from which they result could spoil the outcome 
of nucleosynthesis. We will take care of this, but it will 
not be mentioned in the following. However, on the other 
hand, recent CMB observations allowed to measure the 
baryon density of the universe directly, and some subtle 
strains between this value and the predictions of standard 
nucleosynthesis have shown up |22j. Maybe one reason, 
among others, could be a certain amount of fluctuations 
of the order of the horizon scale at this time. 

For the temperature T k at the time of formation tk , we 
obtain 



H 2 {t k ) 
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=> k B T k = 




(9) 
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Here g(tk) denotes the (effective) number of relativistic 
degrees of freedom at time tk- During the formation 
of the PBHs we are interested in, these degrees com- 
prised electrons, positrons, photons, and three kinds of 
left-handed neutrinos, so we always take g[tk) = 10.75 in 
our calculations. 

The connection between the horizon mass and the scale 
of the horizon is given by the formula 
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=> M H (t k ) = 6.3xl0 12 ^gJ^(^—\ M .(1O) 

At last, for the comoving number density of PBHs of 
mass Mjy(ifc) we get 



n(M H (t k )) = 



(3(tk)p{t k ) fa(t k )H(t k ) 



M H {t k ) V H (tk) 
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It is this formula, together with (0 and ©, which we 
will employ in our considerations. As an input we need 
suitable parametrizations of the primordial fluctuation 
spectrum. A few of them are discussed in the following 
subsection. 



B. Parametrizations of Afj(fc) 

There are several possibilities to parametrize the spec- 
trum of primordial curvature perturbations A|j(fc). Fits 
of observations of LSS and CMB usually employ a power- 
law spectrum with some amount of running: 



A^(fc) 



A~ ( , i)( A 



nQ + h niln {^ 



(12) 



^ n2ln U 



Here kg is some pivot scale, and the parameters no and 
n\ denote the tilt and the running of tilt at that scale, 
respectively: 

rflnA| 
n o -- ,, /, /, x (fco) = ns{k ) , 



"l 



d\n(k/k ) 

dn s 
d\n(k/k ) 



(fc ) = a s (k ) 



(13) 



The simplest models of inflation suggest that the coeffi- 
cients rii scale as powers e l of some slow-roll parameter 
e J2Jj. In this case, for e w 0.1 an expansion up to i = 4 
can be expected to be accurate to 10% for about 16 e- 
foldings around the pivot scale. For a pivot scale of 0.05 
/i/Mpc, this means sensitivity down to horizon masses 
of 400M Q . For general n^ this parametrization can, of 
course, also be seen as a phenomenological, model inde- 
pendent description of the primordial perturbations. 

From the combined analysis of several CMB and Large 
Scale Structure (LSS) obsevations the WMAP team 



claimed to have detected a value of a, 



-0.075 



+0.044 
-0.045 



[12| . The large modulus of this value is quite constrain- 
ing for simple models, and the significance of the obser- 
vation has been doubted by several authors. However, 
recent small scale measurements of the CMB by the Cos- 
mic Background Imager (CBI) yielded further evidence 



for a negative running parameter as — — O.OSTto^g |24j . 
A definitive measurement of as , and possibly also of n<i 
and/or 713, is expected from the Planck satellite in com- 
bination with the full data set of the SDSS and other 
upcoming LSS surveys. 

The spectrum may also be expanded directly in powers 
of the logarithm of the comoving scale 



A\[k) = A 2 R (k ) 




(14) 



This expansion is not suggested by some specific model, 
but it seems to be more appropriate for predictions of 
some more general inflationary models and gives in many 
cases a better reconstruction of the 'real' spectrum than 
(|12f> . cf. I23]. It is clear that cosmological parameter fits 
should be done with a few different parametrizations as 
long as one has no knowledge of the form of the spectrum 
which is realized in nature. Even very accurate data may 
lead to a bad fit if the employed parametrization is not 
suitable. Note that a similar problem occurs when one 
tries to model the time evolution of the equation of state 
of dark energy. 

The Oj can - for fixed fco ~~ be expressed in terms of 
the rii in the following way: 

01 = n -l, a 2 = (n - l) 2 +ni , (15) 

a 3 = (ri - l) 3 + 3(n - l)n\ + n 2 , a 4 = (n - l) 4 
+6(n — l) 2 «i + 4(n - 1)«2 + 3n 2 + n 3 , ... 

However, one should in general not expect that indepen- 
dent fits of the spectra (|T2"|l and (|l4*jl up to a given order 
yield the relations l|15fl . This may happen because the 
number of expansion terms needed to reproduce the real 
spectrum over a certain range of scales is probably dif- 
ferent for parametrizations like (|12f> and l|14|) . 

We now come to a more theory driven and model de- 
pendent parametrization of A^(k). Here it is assumed 
that the fluctuations are generated during an inflationary 
era described by an inflaton potential V((f), the deriva- 
tive of which has a jump at some value </> s corresponding 
to a comoving wave number k s . It has been shown by 
Starobinsky [2|| that the ensuing fluctuation spectrum 
has a universal form around the scale k s , which can be 
calculated exactly and is given by 
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IV. COMPUTATION OF PBH ABUNDANCES 
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FIG. 1: (Color online) Typical form of the BSI spectrum with 
jump parameter p = 0.07. The jump scale k„ is a free parame- 
ter and the large scale normalization may be chosen consistent 
with observational results. 



Here x = k/k s denotes the wave number in units of k s 
and p is the ratio of the left and right handed limits of 
the derivative of the inflaton potential at 4> s . The typical 
form of the spectrum is depicted in Fig. 1. Because of 
the broken scale invariance we will refer to l|16|) as the 
BSI-spectrum. Its asymptotic behavior is given by 



A 2 



A 2 



(k) ™ =£ 
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and approaches the scale invariant Harrison-Zel'dovich 
spectrum on large and small scales, but with different 
amplitudes. If the jump in the spectrum would occur 
at larger scales it could perhaps be observable as some 
enhancement of the small scale structure of galaxy dis- 
tributions, which would be in conflict with actual obser- 
vations. However, the scales we are interested in here are 
too small to be problematic with respect to small scale 
structure. 

The normalization A 2 , may be taken from CMB mea- 
surements, and thus the model has two free parameters, 
p and k s . But even if such a kind of spectrum is realized 
in nature, one would expect additional structure super- 
imposed on it on large scales which depends on the actual 
form of the inflaton potential. 



Now we want to apply the formulae of the last sec- 
tion to see whether PBHs can contribute to the popu- 
lation of seed black holes which then evolve to become 
SMBHs. In (3 the comoving number density of the 
population of z > 6 quasars has been estimated to be 
(6.4 ± 2.4) x 1CT 10 Mpc~ 3 (based on a population of 10 
objects). It is not easy to get accurate constraints for 
the slope of the luminosity function at such an early era, 
but one should expect these objects to be the heaviest 
ones at their time with masses up to about 3 x 10 9 Mq. 
In section two, we have shown that black hole seeds with 
a mass of more than 1300 Mq at z sw 15 are required 
to grow the SMBHs powering these quasars. In order 
to test the hypothesis whether these seeds could be pri- 
mordial, we try to reproduce this density with PBHs. 
Taking into account some inefficiencies in the accretion 
and merging processes involved, we shall study whether 
number densities few orders of magnitude higher than 
the density mentioned above seem to be possible. For 
the masses it should be sufficient to have the peak of 
the distribution at a value of several 100 Mq. A PBH of 
this mass can form a small overdense region of its own 
which is able to grow and take part in merging events 
after the onset of the matter dominated era, according 
to the hierarchical model of structure formation. After 
the nonlinear collapse of the halo the accretion of bary- 
onic matter could start very early and does not have to 
wait until seed black holes have been formed by astro- 
physical processes having their own problems. Another 
mode of early growth could be the "accretion" of energy 
from a surrounding quintessence field |35| |27| or the swal- 
lowing of interacting dark matter. Thus, it seems to be 
possible to have seed black holes of the desired mass at 
z sw 15 from PBHs which are slightly lighter, and we will 
be content if the latter have masses of several 100 Mq, 
as mentioned above. 

We start with a very simple spectrum of primordial 
curvature fluctuations which, nevertheless, is up to now 
consistent with all observations: a scale free power law 



A 2 R (k) = 2.95 x 10" 
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(18) 



A fit of several recent CMB and LSS observations at the 



•0.019 



pivot scale fc = 0.05 Mpc" 1 yields^ A = 0.631; 

and ns — 0.966t2!o20 U3- Putting these values in the 
formulae of the preceding section gives - as could be ex- 
pected - a rather disappointing result: The expected 
number density of black holes is absolutely negligible. 
This is not surprising given the known similar results of 



attempts to establish a significant contribution of dark 
matter from PBHs [ll|. So we arrive at our first result: 
// the spectrum of primordial curvature fluctuations from 
scales relevant for CMB and LSS up to k m 10 5 Mpc~ 
is given by a scale free power law with parameters con- 
sistent with CMB and LSS fits, then the seed black holes 
for future SMBHs cannot be primordial. 

It does not make much sense to ask which constant 
parameter n$ would be needed in order to produce an 
appreciable amount of PBHs with the desired mass for 
the following reason: The number density is in this case a 
sharply decreasing function of the PBH mass, and so one 
runs very quickly into an overproduction of small mass 
PBHs. These would overdose the universe for n s > 1.35 
|28|. If one cuts off all fluctuations on scales smaller 
than the ones needed to produce PBHs of about 1000 M Q , 
a value of n s ss 2.1 would be required to produce the 
desired number of black hole seeds. 

It is clear that an extrapolation of a scale free power 
law over such a large range is highly questionable. Note 
that this remark applies in particular to the attempts to 
compute dark matter abundances with such a spectrum, 
because the PBH masses which seem to be preferred in 
such studies are of the order 10 15 — 10 20 g and are there- 
fore related to scales much smaller than the ones consid- 
ered in the present work. 

Let us now take some more sophisticated parametriza- 
tions of primordial curvature fluctuations. We start with 
the "running tilt" -spectrum (|12fl with a pivot scale of 
fco = 0.05 Mpc -1 and values of the normalization con- 
stant as well as no and n\ consistent with cosmological 
fits. After some trial and error we see that it is possible 
to find values for the remaining two parameters which 
describe a fluctuation spectrum that is able to produce 
a reasonable number of PBHs of the desired mass. The 
result can be seen in Fig. 2, where we have plotted the 
expected comoving number density of PBHs per e- folding 
in M, which is more intuitive than the spectral number 
density. Although the adopted parameter values do not 
seem to be completely unphysical, one should be careful 
to give them (or the amplitude and scale of the bump 
in the fluctuation spectrum they lead to) too much rel- 
evance. We only want to show that such a spectrum is 
still possible. It is important to note that it should be 
feasible for future observations to obtain a quite accurate 
value for n\ and possibly more or less strong constraints 
on U2- The contributions of the terms of order five and 
larger in the expansion of the fluctuation spectrum are 
irrelevant for the scales wc arc interested in, provided, of 
course, that the parameters are of the order of unity at 
most. Since the resulting PBH density is quite sensitive 
to U2, it might be possible to falsify the idea of having 
PBHs as SMBH seeds. 

One might think that the message of the previous para- 
graph is by no means surprising, and that for any kind 
of parametrization of primordial fluctuations some clever 
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FIG. 2: (Color online) Expected comoving number density 
of PBHs per e-folding in M for the spectrum of primordial 
curvature fluctuations 11211 with adopted parameters no = 
0.98, ni = -0.05, n 2 = 0.123, n 3 = -0.022, and all other m 
equal to zero. 



choice of parameters could be found which does the de- 
sired job. However, this is not the case as can be seen in 
the example of the parametrization (|14|) . Because we are 
not aware of any published fit of its parameters, we will 
take the values a\ — —0.03 and 02 = —0.05 in the follow- 
ing, which are suggested by the relations (|15fl . Playing 
the same game as above, it turns out to be impossible 
to find satisfying higher order parameters to produce an 
appreciable amount of PBHs. More precisely, we need 
more parameters than can ever be sensibly fitted, and the 
values of these parameters are unpleasingly large. This 
means that if the parametrization (|14fl , together with the 
values of a\ and 02 mentioned above, should ever lead to 
a good fit of the primordial fluctuations, then PBHs are 
strongly disfavored to be the seeds of SMBHs. 

The reason for this is that we have assumed the first 
coefficients a\ and 02 to be negative. Thus less fluctua- 
tions are expected on small scales than on large scales. 
To produce PBHs we have to achieve a turnover for which 
quite large positive values of the following parameters are 
needed. But, unfortunately, these overproduce PBHs on 
smaller scales and so further negative parameters have to 
be introduced. Although the parametrization (|12f) is also 
plagued by this problem, it turns out to be not so dra- 
matic as in the present case because there the expansion 
is performed in the exponent. 

The conclusions of the preceding paragraph relied cru- 
cially on the fact that the parameters a\ and 02 were 
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as a more convenient parameter. The maximum of the 
distribution shows up at a slightly lower scale than M s , 
which is not surprising because the bump in the fluctua- 
tion spectrum also appears at scales smaller than k s (see 
Fig. f). 

Although it seems to be possible to obtain the de- 
sired number density of PBHs with the BSI-spectrum, it 
should be noted that in Fig. 3 we have chosen for the large 
scale normalization of the fluctuation spectrum a value 
which is situated at the upper end of its confidence re- 
gion and, more important, the adopted value p — 0.0004 
leads to fluctuations of the order 1% at the bump, which 
itself is quite narrow. It does not make sense to consider 
smaller values of p in our calculations, because then lin- 
ear theory would break down and we could not trust our 
results. Furtheron, it is hard to accept such fluctuations 
from a realistic point of view. 

The BSFspectrum could be further constrained by an- 
alyzing specific inflaton potentials, which would lead to 
further structure in addition to the universal part of the 
spectrum, in particular on large scales, where it could be 
better anchored to observations. 



FIG. 3: (Color online) Expected comoving number density 
of PBHs per e-folding in M for the BSI-spectrum of pri- 
mordial curvature fluctuations 116H with adopted parameters 
p — 0.0004, M 3 = 2OOOA/0, and large scale normalizationx 
Ag = 3xl0 -9 . 



assumed to be negative. Hovewer, present observations 
are not yet able to exclude positive values for them com- 
pletely, and in this case it would be easier to extend a 
good fit of (|T4"Jl to a spectrum suitble to produce a con- 
siderable number of PBHs, but, of course, a good fit of 
(|12fl with a positive value of as would still be more fa- 
vorable for our intentions. 

Now we come to the parametrization of the primordial 
fluctuations according to the BSI-spectrum 1|16[) . Here 
we have only two parameters, p and fc s , at our disposal, 
and it seems to be much more difficult to produce PBHs 
out of this spectrum. But in practice things turn out 
not to be that worse because of the following reason: 
what we need is a spectrum which is compatible with 
observations on large scales, shows a large bump at a 
scale suited for producing PBHs of the desired mass, and 
has less fluctuations on smaller scales in order to avoid 
an overproduction of small PBHs. As can be seen in Fig. 
1, the BSI-spectrum has such a form, and the parameters 
p and k s are exactly the screws for the position and the 
amplitude of the mentioned bump. It should be noted 
further that, contrary to the preceding parametrizations, 
the spectrum is fixed to its known large scale behavior by 
only one parameter, namely the overall normalization. 

A result of our juggling with the parameters can be 
seen in Fig. 3. Instead of the scale k s we have introduced 
via eq. (|10|) the PBH mass M s connected with this scale 



V. CONCLUSIONS 

We have studied whether the black hole seeds needed 
to grow SMBHs could be primordial ones formed within 
the very first second of cosmic evolution. In order to do 
this the expected number densities for PBHs of the rel- 
evant mass (« 1000 M Q ) have been calculated for a few 
parametrizations of the spectrum of primordial curvature 
fluctuations normalized to their known large scale val- 
ues. As expected, it is impossible to produce the required 
black holes with a scale free spectrum of the form i|18|) 
without a huge overproduction of smaller PBHs which 
eventually closes the universe. 

What is required, is a fluctuation spectrum which is 
consistent with observations on large scales, shows a large 
bump with amplitudes up to several parts of a percent 
on scales of about O.lpc -1 , and decreases quickly on 
smaller scales in order to avoid an overproduction of 
PBHs. This became clear in the considerations of three 
different parametrizations which seem to be capable of 
achieving the mentioned demands. Two of them were 
expansions in powers of In A:. In the first one the expan- 
sion was done for the exponent of a power law, and in 
the second one for the fluctuation spectrum itself. Two 
parameters, the overall normalization and the power law 
slope at some large pivot scale, can be fixed quite ac- 
curately with present observations. A third one is only 
vaguely constrained today. It is, of course, possible to 
reproduce the desired form of the spectrum with an ar- 
bitrary number of parameters, but if one restricts the 
considerations to parametrizations which employ only a 
small number of parameters in addition to the ones fixed 
by observations and postulates further that their modu- 
lus should not be larger than about unity, as is suggested 
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by inflationary models, than the scenario of having PBHs 
as SMBH seeds could be strongly disfavored in the near 
future. This will be the case for the second parametriza- 
tion if the evidence of a negative value of the running 
parameter as should be confirmed. 

The third parametrization, the BSI-spectrum, on the 
other hand, employed only two adjustable parameters 
and one parameter, the overall normalization, to anchor 
it to its large scale value. As a universal spectrum it can- 
not be further constrained by observations, but it is clear 
that there must be some superimposed structure on top 
of it, which arises from the specific form of the inflaton 
potential. This might be delineated by future observa- 
tions, and so the position of the characteristic scale of 
the spectrum could be constrained, which then perhaps 
leads us to abandon that model. 

It seems to be quite difficult to construct viable phys- 
ical models which produce the desired jump in the spec- 
trum of primordial fluctuations. However, we have al- 
ready seen that the BSI model does the job. Other sce- 
narios are, e.g., the formation of domain walls during 
inflation |29| , one or more phase transitions during infla- 
tion |30| , or some events of resonant particle production 
in the inflationary era [31|. Most of these mechanisms 
have been invented as a mean to produce features in the 
primordial fluctuation spectra which are detectable by, 
e.g., CMB, LSS or Ly-a observations. However, there 
seems to be no reason why they should not be applicable 
on smaller scales, which are relevant for the considera- 
tions of the present work. 



It should be noted that the required scale of the bump 
in the fluctuation spectra is near the scale relevant for 
fluctuations triggered by the quark-hadron phase transi- 
tion [33. However, the PBHs resulting from this fluc- 
tuations are expected to be lighter than « 1M©, and 
so a larger amount of early growth would be required, 
perhaps by the absorption of a quintessence field 
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Future observations and experiments promise to be 
able to decide whether PHBs are serious candidates for 
SMBH seeds. High energy physics may shed some light 
on the specific mechanism of inflation, in particular the 
on scale of reheating. CMB, LSS, and weak lensing ob- 
servations are expected to proceed in reconstructing the 
spectrum of primordial curvature fluctuations to smaller 
scales. The epoch of galaxy formation and reionization 
will be reached with the next generation of telescopes, the 
space based gravity wave detector LISA will clarify the 
abundance and importance of mergers of SMBHs, and 
further supernovae observations should explore the evo- 
lution of dark energy, which could lead to a reassessment 
of the formation times of the early generation of quasars. 

Thus, primordial black holes, which have been a matter 
of debate for decades, may, at the end, show up at places 
where they had hardly been expected! 
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